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Frequency Analysis of Stationary, 
Non-stationary and Transient Signals 


Introduction 


Frequency analysis is a very im- 
portant technique in investigations 
on many types of physical systems 
such as mechanical, acoustical and 
electrical systems. One of the main 
reasons for this is that many such 
systems have approximately linear 


~ properties, which means that if they 


are excited at a certain frequency, 
then the response occurs at the 
same frequency, though with ampli- 
tude and phase modified according 
to the transfer function. As a result, 
the response at this frequency is in- 
dependent of what is happening at 
other frequencies, and thus a de- 
composition of the signals into their 
components at various frequencies 
represents an appreciable simplifica- 
tion. 


Moreover, the system response 
(transfer function) can most easily 
be expressed as a function of fre- 
quency since the combination of 
cascaded systems then _ involves 
only multiplication of their  fre- 
quency characteristics (simplified to 
addition on taking logarithms) com- 
pared with convolution of the equiv- 
alent characteristics in the time do- 
main. 


The basis of frequency analysis 
are the Fourier integral pair: 


F(o) = = i f(tje—iot dt (1) 
and 
f(t) = i. F(w) eiwt dw (2) 


by R. B. Randall, B. Tech., B.A. 


An interpretation of the first equa- 
tion is that for a time function f(t), 
the complex frequency component 
F(w) is obtained by correlation of f(t) 
over its entire length with the rotat- 
ing unit vector e—/#t (with angular 
velocity —w), while the second 
states that f(t) can be regenerated 
by summing a (possibly _ infinite) 
number of rotating vectors which at 
time zero are equal to F(w) and 
which rotate with angular velocity 
Ww. 


In many cases, phase _ relation- 
ships are of secondary importance, 
and it is only necessary to find the 
amplitude or energy distribution 


with frequency. In such cases, 
Parseval’s relationship is directly 
relevant: 

f [Alt]? dt= S| Flw)|? deo (3) 


which states that the total energy 
can be obtained either by integration 
over all frequency or over all time, 
and also that the energy density at 
angular frequency w is given by 
| F(w)|2 where F(w) is the complex 
component defined in (1). 


The mathematical expressions 
(1), (2) and (3) assume that the 
ranges of integration are infinitely 
long and also that the analysis band- 
width dw is infinitely narrow, which 
of course cannot occur in practice. 
Moreover, no mention has been 
made of the mathematical condi- 
tions for the existence of the inte- 


grals, but these are also connected 
with the infinite ranges of integra- 
tion which do not exist in practice. 


The purpose of this Application 
Note is to explain how practical vat- 
ues can be obtained for | F(w)| (gen- 
erally expressed on a logarithmic or 
dB amplitude scale) using different 
types of analyzer (having a finite res- 
olution bandwidth) and for different 
types of signal. The latter include 
both stationary (periodic, quasi-peri- 
odic and random) signals and non- 
Stationary (continuous and _ tran- 
sient) signals. 


The discussion includes such prac- 
tical questions as choice of band- 
width, averaging time and analysis 
speed, as well as the use of special 
weighting functions or ‘time win- 
dows” to select out a particular part 
of a continuous record and/or to 
avoid problems with the discontinul- 
ties at the ends of a signal seg- 
ment. 


A discussion is also included of 


- the various possibilities for the anal- 


ysis of transients (including win- 
dowed segments). These can be 
joined into a loop and analyzed as 
periodic signals or alternatively the 
energy content at each frequency 
can be measured by passing the 
transient once through the filter for 
each filter position. 


The treatment is of a Summary na- 
ture and refers back to the original 
sources for detailed derivations. 


Analysis Methods 


There is a choice between the fol- 
lowing types of analyzer: 


Fixed Filter 

This type of instrument has a 
number of fixed contiguous filters 
which together cover the whole of 
the frequency range of interest. A 
detector circuit successively meas- 
ures the averaged output from 


each filter in turn, and the results. 


would normally be recorded on a lev- 
el recorder. This type of analyzer is 
generally only used for relatively 
broadband analysis, such as 1/3 oc- 
tave, since it would otherwise re- 
quire too many filters. The B&K 
Spectrometers Types 2113 and 
2114 are typical. 


Sweeping Filter 

This type has a filter with a vari- 
able centre frequency which is 
Swept over the frequency range of 
interest. There are two main types, 
with constant absolute bandwidth 
and constant relative (percent) band- 
width respectively. The B&K Ana- 


lyzer Type 2010 is typical of the first’ 


type, and operates according to the 
heterodyne principle. It should be 
kept in mind, however, that the 
bandwidth can increase stepwise 
(between 3,16 and 1000 Hz) with in- 
creasing frequency, thus approxi- 
mating a constant proportional band- 
width filter. 


The Analyzers Types 2120 and 
2121 have constant proportional 
bandwidth, but this can be selected 
over a wide range between 1% and 
23% (1/3 octave). The older Ana- 
lyzer Type 2107 is also of this type. 


High Speed Analysis 

When a wide frequency range is 
to be swept with a relatively narrow 
bandwidth, then normal swept fre- 
quency analysis can take a _ very 
long time. In this event a consider- 
able time saving can be made by re- 
cording the signal on a Digital Event 
Recorder (Type 7502) and playing it 
back at a considerably higher 
speed. As will be shown, this per- 
mits a correspondingly faster analy- 
sis. The analysis time can typically 
be reduced from over an hour down 
to a few minutes. 


Real-Time Analysis 

This means basically that the anal- 
ysis is so fast that it takes less time 
than that required to collect the 
data record analyzed. A number of 
different methods are used _ to 
achieve this, even up to frequencies 
as high as 10 to 20 kHz. 


Analogue Parallel Analysis is the 
principle employed in the B&K An- 
alyzer Type 3347 (1/3 octave band- 
width). This not only has a number 
of separate filters (like the spec- 
trometers) but also a separate detec- 
tor for each filter, so that the output 
levels from all filters are displayed 
simultaneously on a screen. Once 
again, the number of separate fil- 
ters and detectors sets a limit as to 
how narrow band the analysis can 
be, and the principle is normally not 
used for bandwidths less than1/3 
octave. 


The Time Compression _ tech- 
nique is used in the B&K Analyzer 
Type 3348 which has 400 narrow 
(constant) bandwidth channels. The 


basic principle is similar to that de- 
scribed under ‘‘High Speed Analy- 
sis’’ except that recording in the dig- 
ital memory goes on during the 
same time as playback and hetero- 
dyne analysis. The frequency trans- 
formation ratio is also much higher 
which means that real-time analy- 
sis is possible up to 10kHz. Con- 
nected to the analyzer unit is an av- 
eraging unit which can carry out av- 
eraging over a selectable number of 
Statistically independent spectra. 


Digital Analysis should be men- 
tioned to complete the picture. Two 
principles are used, viz. digital fil- 
ters and FFT (Fast Fourier Trans- 
form) techniques, both of which can 
be realized either as software or 
hardware. The latter would _ nor- 
mally be required for real-time anal- 
ysis, however, and hardware instru- 
ments have so far been relatively ex- 
pensive. | 


Digital filters are similar to ana- 
logue filters in that the math- 
ematical description of the latter 
is converted to a numerical algo- 
rithm with which the input values 
(from an A/D converter) are pro- 
cessed. The RMS values of the out- 
put series are also calculated nu- 
merically. 


FFT is a very efficient algorithm 
which calculates the DFT (Discrete 
Fourier Transform), the discrete ver- 
sions of equations (1) and (2). It 
makes possible a very rapid transfor- 


mation (in either direction) between 


time and frequency domains, with 
retention of phase information if so 
desired. 


Choice of Analysis Parameters — Analysis Speed 


This section is mainly concerned 
with swept frequency narrow band 
analysis of stationary _ signals, 
where considerable time can be 
saved by analyzing as efficiently as 
possible. 


The parameters which influence 
analysis speed are the following: 


Bandwidth 

This in fact has a double effect, 
since it is obvious that the analysis 
time is proportional to the number 
of bandwidths in the frequency 
range to be covered. In addition — 
as we shall see — the bandwidth it- 
self normally also has an effect on 
the time to be taken in sweeping 


over a bandwidth. Thus it is import- 
ant to use the largest possible band- 
width which still gives fully useable 
results. 


The first choice to be made is be- 
tween constant absolute and con- 
stant relative bandwidth. It is not 
possible to give definite rules for 


this, but the following discussion 
may be helpful: 


It is often claimed that “narrow 
band analysis” is virtually synony- 
mous with constant bandwidth anal- 
ysis, but this is not necessarily the 
case. As an example, the Analyzer 
Type 2120 has bandwidths down to 
1% and this will often give suffi- 
cient resolution. In fact, many physi- 
cal systems behave in principle like 
constant relative bandwidth filters 
(constant Q-factor being equivalent 
to constant percentage bandwidth) 
and thus where the excitation is 
fairly broadband, it will be most effi- 
cient to analyze the response with 
constant percentage filters. 


A further advantage of constant 
relative bandwidth is that it gives 
constant resolution on a logarithmic 


frequency scale and can therefore 


be used over a wide frequency 
range. Other grounds for using a log- 
arithmic frequency scale (though 
“not necessarily constant percentage 
bandwidth) are: 


(a) Small speed changes in, say, a 
machine only give a lateral displace- 
ment of the spectrum. 


(b) Certain relationships can most 
easily be seen on log-log scales. 


Where Q-factors exceed about 50 
it may in any case be necessary to 
employ constant absolute band- 
width purely to obtain less than 1% 
bandwidth in a part of the fre- 
quency range. (Note, however, that 
the already-mentioned stepwise ap- 
proximation of constant percentage 
bandwidth offered by the 2010 
does not suffer from this’ limita- 
tion). 


Constant absolute bandwidth (in 
particular with a linear frequency 
scale) is normally most beneficial 
when the signal consists primarily 
of a number of harmonically related 
components, since these are then 
equally separated and_ resolved 


(Fig. 1). 


Integration Time 
For gaussian random signals the 
accuracy of the result is given by 
the following equation (Ref. 1). 
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Fig.1. Vibration spectrum with a large number of harmonically related components 
where 
€ = standard deviation of the error lyzer averaging time (T,) or the 
B = filter bandwidth (Hz) length of the signal sample 
Te = equivalent integration time de- (Ts) according to the following 


termined either by the ana- 


Linear 
Integration 


table: 


Exponential (RC) 
Integration 


Ta 
Ts 
0,75T, (= 0,75Ts) 


Table 1. 


As an illustration, 


| 


€ ~ 1,5dB for BTe = 10 
and 
€ ~0,6dB for BT¢ = 50 


For deterministic signals (periodic 
Or quasi-periodic, made up entirely 
of sinusoidal components) the only 
requirement of the averaging time 
(for reasonable accuracy) is that it 
Spans over approximately 3 periods 
of the lowest frequency present (or 
7 periods in the case of half-wave 
rectification) in order to reduce the 
ripple to- less than +1/4dB. It 
should normally be possible, how- 
ever, to increase the averaging time 
to a half of the actual dwell time 
per bandwidth (Tp) without having 
to reduce analysis speed. If higher 
accuracy than about 1dB is de- 
sired, then the averaging time 
should be related to the filter re- 
sponse time (see later) by making 
BT pa 210, Say. 


For signals on a loop (tape loop, 
or digital event recorder) the averag- 
ing time must extend over a num- 
ber of loop lengths (at least 3) in or- 
der to achieve a ripple-free result. A 
lower averaging time can however 
be used (see Table l) if the resulting 
fluctuations do not confuse the re- 
sults (e.g. by neglecting the splice 
noise peaks which occur with a 
tape loop). 


Averaging time in connection 
with impulse analysis is discussed 
in detail later. 


Filter Response Time 

When analyzing continuous deter- 
ministic signals, the sweep speed 
will normally be determined by the 
filter response time (Tp) which is 
given approximately by 


1 
Tr ~B seconds (5) 
Fig.2 shows a typical filter re- 


sponse (to a sinusoidal input) for a 
1/3-octave filter. 


Filter Input Signal 
0 
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Fig.2. Typical filter response from a 1/3 octave filter 


Sweep Speed with DC-Recording 

The newer and more flexible in- 
struments (2010, 2120, 2114) 
have a DC output with 60dB dy- 
namic -range, giving the possibility 
of using the internal detector even 
for recording the analysis on a level 
recorder (which is then used in DC 
mode). As will be shown, the record- 
er does not then have any limiting 
effect on the sweep speed. 


Defining the filter dwell time (Tp) 
_ by the following equation: 

B 
Sweep Speed = Tp (6) 


Then Tp is determined as the 
greater of either 


Tp 2n1 Ta (7) 
(with respect to averaging time) 
or 

Tp 2 no Tr (8) 


(with respect to filter response) 


* Defined in Fig.3. 


integration integration integration 
1,2 dB 0,5B 0,3B 
n> = 10,0. =1dB 


Table 2. 


Equation (7) normally governs in 
the case of random signals and sig- 
nals on a loop, while (8) governs in 
the case of continuous deterministic 
signals (accuracy ~ 1 dB). 


The values of n, and n», are de- 
pendent on the desired accuracy of 
measurement, but typical values 
are given in Table 2 (valid for spec- 
trum peaks). 


Further information — including 
data on valleys in the spectrum is 
given in Ref.2. 


It can be seen from equations (4), | 


(5), (7) and (8) that Tp is always gov- 
erned by an expression of the type 
B Tp 2K. e.g. for a random signal 
with BT- = 10 (Equ.4) Te = Tay 
(Table 1) and n, = 1-— (Equ.7, 
Table 2), BT, must be 2 10. 


Alternatively, for a deterministic 
signal with BTy ~1 (Equ.5) and n> 
4 (Equ.8, Table 2), BTp must be 
4. 
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Fig.3. Error and delay in writing out peaks 
in a spectrum 


Fig.4 gives lines corresponding to 
the most commonly used values of 
K for B&K’s standard bandwidths 
(both constant and constant percent- 
age) so that it is possible to select 
the correct value of Tp for each fre- 
quency range (Ref.3). The reference 
also includes tables indicating the 
equivalent of the bandwidth in mil- 
limeters on the recorder paper, so 
that the paper speed can be deter- 
mined directly by dividing this with 
Tp- 


The fact that constant accuracy re- 
quires a constant BTp product is 
the explanation of why frequency 
transformation upwards allows the 
analysis to be carried out faster. 
Both the total analysis range and the 
analysis bandwidth are increased by 
the same factor (for unchanged reso- 
lution) but the increased bandwidth 


permits a corresponding decrease of 


Tp. Since this is the sweep time 
per bandwidth, the total analysis 
time is also reduced in the same 
proportion. 


AC Recording 

In some cases (e.g. no DC output 
from the analyzer, desired dynamic 
range >60dB, desired integration 
time <0O,1s) the AC output from 
the analyzer is taken directly into 
the level recorder which then per- 
forms the averaging. The effective 
averaging time depends on the writ- 
ing speed settings, as given in 


Ref.4 (conservative values also 
given in Fig.4). However, in this 
case it will often be this writing 


speed which limits the analysis 
speed, depending on its ability to 
write out the spectrum. The worst 
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For each bandwidth and frequency range, determine Filter Dwell Time, Averaging Time and Sweep Speed as follows: 


1. Random Signal: Dwell time Tp is always determined by averaging time T,. 
From the graph read off Tp ( = Ta) from the appropriate line for BT, = 10 (horizontal lines for constant bandwidth, sloping lines for 
constant proportional bandwidth). This value will correspond to + 1 + dB RMS error. For a higher BT, product and consequent re- 
duced error, increase the value of T, proportionally. 
For DC recording calculate sweep speed S = tT): 
For AC recording, T, is determined by recorder writing speed W and this can be read directly from the right-hand scale. Sweep speed can 


then be calculated from S = (applicable to 50 dB potentiometer and 100 mm acl 


100 


2. Deterministic Signal (periodic or quasi-periodic): Read Tp based on filter response time from appropriate line for BTp = 4 (or 2 where 
bandwidth < 1%). Read also minimum averaging time T, based on + 1/4 dB ripple from line fT, = 3 (independent of bandwidth). 


For DC recording calculate sweep speed based on filter response as S = = 
D 


For AC recording read off writing speed W corresponding to T, and calculate sweep speed based on recorder response as S = ~~. The 
governing sweep speed is the lesser of this and the value as calculated for DC recording. 


In all cases where sweep speed is governed by filter response time, increase Ta to the highest value which does not affect sweep speed. 


Note (1) Averaging times will normally be constant over at least a half decade, whereas filter response times change automatically with frequency. 


Note (2) S will be in Hz/s for B in Hz, but will be directly in mm/s for B expressed as equivalent recorder paper length in mm. 


272264 


Fig.4. Basic analysis parameters 


case is represented by recording the 
filter characteristic, and the steep- 
est slope of filter flank likely to be 
encountered is 50OdB/bandwidth 
(shape factor 3,4). Simple geomet- 
ric considerations (Fig.5) show that 
the sweep speed will then also be 
limited by: 
- BW 


SH Te) (9) 


S = sweep speed 


and W 


writing speed (mm/s) 


S will have the units Hz/s when 
B is expressed in Hz, and mm/s for 
B expressed in mm. 


When Tp is determined by Ta 
(random signals, looped _— signals) 
equ. (9) will always govern, but 


when Tp is determined by Tr (con- 
tinuous deterministic signals) it has 
to be checked whether equ. (9) or 
(6) governs. 


In DC-recording mode it is always 
possible to use 1000mm/s writing 
speed (2000 mm/s for Level Record- 
er Type 2307) which when substi- 
tuted in equ. (9) results in: 


S<10B 


Furthermore, for DC recording the 
lowest possible Tp is equal to the 
minimum existing averaging time 
T, (0,18) which when substituted 
into equ. (6) also results in 


S<10B 


which shows that equ. (9) will 
never be a limiting factor in DC-re- 
cording. 


High Speed Analysis 


As mentioned previously it is pos- 
sible by recording a signal at one 
speed and playing it back at a 
higher speed, to reduce the analysis 
time accordingly. This applies 
equally for tape recorders, digital 
event recorders, and time compres- 
sion analyzers. The latter two give 
much more scope, however, since 
digital techniques allow much grea- 
ter frequency transformations with- 
out the mechanical problems which 
limit tape recorders. 


In order to reproduce these flanks, 
paper speed limited by 


P< an mm/s 


Writing Speed 
W (mm/s) 


<—__—___————Paper Speed P (mm/s) 


740305 


Fig.5. Limitation of recorder writing speed in AC recording 
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Fig.6. (a) Normal swept analysis 
(b) High speéd analysis with rectangular window (10 K) 
(c) High speed analysis with gaussian window 


This means, however, that it is 
necessary to take a segment out of 
the continuous’ signal which is 
joined into a loop and analyzed as a 
periodic signal (with period equal to 
the loop length). This has two main 
effects: 


(a) The spectrum of a periodic sig- 
nal consists of spectral lines with a 
spacing equal to 1/T (where T = the 
period). The heights of the spectral 
lines are however proportional to 
the corresponding continuous spec- 
tral density function, which can 
thus be considered to be “sampled” 
at the discrete frequencies. It can 
be seen, however, that both the 
minimum valid frequency and the 
minimum valid bandwidth are lim- 
ited by the line spacing 1/T. On 
the other hand the upper frequency 
is limited by the antialiasing filters 
before the A/D converters. 


(b) A random discontinuity occurs 
where the ends of the segment are 
joined into a loop. This can give rise 
to spurious components not present 
in the original signal, but when the 
signal segment is relatively long, as 
with a 7502 with 10K memory, 
then the effect of the discontinuity 
can often be disregarded. 


This section is primarily con- 
cerned with high speed analysis us- 
ing only a 7502 for frequency trans- 
formation. The next section dis- 
cusses how ‘time windows’ can be 
used to minimize the effect of the 
discontinuity when this is neces- 
sary. 


Provided the analyzer bandwidth 
is greater than the Fourier line spac- 
ing 1/T (which it should be to make 
the spectrum continuous again) 
then the only effect of the discontin- 
uity is a broadening of the base of 
the overall filter characteristic, re- 
sulting in a reduction of selectivity 
(Fig.6). 


Information about the resulting se- 
lectivity is given in Ref.5 for both 
constant bandwidth (as shape fac- 
tor) and constant proportional band- 
width analysis (as octave selectiv- 
ity). Figs.7 and 8 give the final re- 
sults in graphic form. It can be seen 
from Fig.6 however, that the re- 
duced selectivity will seldom be a 
significant limitation in analysis of 
stationary signals. 
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Fig.7. 40 dB Shape Factor for rectangular window length T and ideal filter of bandwidth B 
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Fig.8. Octave selectivity for rectangular window length T and ideal filter with constant percent- 


age bandwidth 


The procedure is quite simple. 
The signal is recorded on the 7502 
at the lowest sampling rate which 
still retains the highest frequency of 
interest (low pass filter cutoff = sam- 


pling frequency/4). The signal is 
then played back at the highest sam- 
pling rate permitted by the 7502 or 
the analysis instruments. In some 
cases the frequency transformation 


will be limited to even decades so 
that the final equivalent bandwidth 
can be made equal to the original. 


Two commonly used possibilities 
are as follows: 


Constant bandwidth (Analyzer 
2010) (Fig.9) 

In this case the maximum play- 
back speed of the 7502 (500kS/s) 
can be used, since the 2010 can 
analyze the maximum resulting fre- 
quency of 125kHz. With a 10K 
memory, the loop length (T) is then 
20ms, so that the minimum avail- 
able averaging time (100ms) can 
be used to smooth out the fluctua- 
tions. The minimum valid band- 
width (1/T) = 50Hz which means 
that there is a choice between the 
following possibilities: 


Filter bandwitdth (Hz) 


100 316 1000 


1250 400 125 
245 74 245 
2 63 20 


Equiv. No. of lines 
Analysis time (s) | 
BTe product (Equ.4) 


Table 3. 


Constant proportional bandwidth 
(Analyzer 2120) (Fig.10) 
In this case the playback speed of 


the 7502 is limited to 100kS/s by | 


the maximum analysis frequency of 


Bandwidth (% 


millimeter equivalent (min.)} 0,5 1,5 5 12 
3 10 10 
250 83,3 25 25 
6 20 46 


paper speed (mm/s) 1 


Analysis time/decade (s) 
BTE product (min.) 2 


Table 4. 


Comparison of Normal, 


It is worth comparing high speed 
analysis with both normal and real- 
time analysis (using the time com- 
pression principle) since it lies some- 
where in between with respect to 
cost and speed. 


(1) High Speed vs. Normal 
Analysis 


Advantages 

Higher speed — narrower band- 
width possible with constant band- 
width analysis — field recording 


time can be much shorter than the 
normal analysis time — short signal 
samples can be analyzed — possi- 
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Fig.9. Typical instrument set-up for high speed analysis with 2010 (constant bandwidth) 


Fig.10. Typical instrument set-up for high speed analysis with 2120 (constant proportional 
bandwidth) 


the 2120 (20kHz). Even so, a con- 
siderable benefit is_ still derived 
from always analyzing in the fre- 
quency range 2 to 20 kHz with fixed 
analysis parameters, with respect to 
both speed and simplicity. At a sam- 
pling rate of 100kS/s (and 10K 
memory) the loop length is 100 ms, 
requiring an averaging time Ta = 
0,3s to smooth out the fluctua- 
tions. 


The figures in Table 4 apply to an- 
alysis made on the standard record- 
ing paper OP 1130. 


For example, to analyze over the 
3 decades 2 Hz— 2 kHz, the record- 
ing sample rate is set successively 
to 0,1, 1 and 10kS/s, while the 
analysis is always performed be- 
tween 2 and 20kHz with playback 
speed 100 kS/s. 


High Speed & Time Compression Analysis 


ble to frequency transform up into 
the frequency range of the analyzer. 


Disadvantages 

One more instrument required — 
both minimum resolution and accu- 
racy limited by sample length. It 
should however be kept in mind 
that a high speed analysis system 
includes a normal analysis system. 


(2) High Speed vs. Time Com- 
pression Analysis 

Advantages 

More flexible (e.g. constant and 
constant percentage bandwidth, lin- 


ear and_ logarithmic frequency 
scales, variable resolution) — the 
various instruments can be used sep- 
arately — better resolution (up to 
2500 lines) — lower cost (espe- 
cially when the normal analysis sys- 
tem already exists). 


Disadvantages 

Not real-time — only means of in- 
creasing the BT¢ product is by in- 
creasing bandwidth and _ reducing 
resolution (however, a BI¢ product 
of 10 is achieved with a 250 line 
analysis) — not so simple to oper- 
ate. 


Effect of Time-weighting Functions 


Taking a segment out of a continu- 
ous signal can be interpreted as 
multiplying the signal by a weight- 
ing function (or viewing it through a 
“time window”). This weighting 
function would normally be rectan- 
gular (when the rest is simply cut 
off) but other rounded weighting 
functions can be used to advantage 
to nullify the discontinuity (Fig.11). 
As shown in Ref.5, the effect can 
also be interpreted as a filtering in 
the frequency domain with a filter 
characteristic equal. to the fre- 
quency spectrum of the weighting 
function itself (Fig.12). It should be 
pointed out, however, that all func- 
tions other than a rectangle give a 
different weighting to different sec- 
tions of the segment, as well as in- 
creasing the minimum valid band- 
width. (Ref.5). 


Fig.13 compares the spectra of 
the following window _ functions 
with that of a rectangle with the 
same total length T (for the gauss 
window interpreted as 3 times the 
half amplitude length as for the 
Gauss Impulse Multiplier Type 
5623). The various time windows 
have the following properties. 


Rectangular 

Minimum ‘possible bandwidth 
(1/T) — highest sidelobe at -13dB 
— sidelobes fall off at 20dB/ 
decade. 


Gaussian function 
Bandwidth ~2/T — no sidelobes 
and very steep characteristic. 


Hanning 

Bandwidth ~1,5/T —_. highest 
sidelobe at —33dB — sidelobes fall 
off at 60 dB/decade. 


Hamming 

Equivalent to a Hanning window 
on a small rectangular base. This 
cancels the highest sidelobe of the 
Hanning function (highest now at 
—43 dB) but means that the remain- 
ing sidelobes. only fall off at 
20 dB/decade. Bandwidth ~ 1,4/T. 


In many cases the Gaussian win- 
dow is the best possible since it 
gives the least disturbance in combi- 
nation with another filter (e.g. 
2010) while the other two are often 


Discontinuity 


PV VW 


Gaussian Impulse 


Trigger 
| 


Delay Time 


Pi signal multiplied by gaussian impulse 


yo 


Discontinuity nullified 


Fig.11. Effect on the junction discontinuity of a gaussian window (time domain) 


Time Domain Frequency. Domain Result from Sweeping 
(Power Spectrum) Ideal Filter 


at 


—co <g— — 3 co 


= Qnf_tt 
p(t) = cos (27 f, t+@,) Pf) 


act 
ty + fo 
dB 


: T als 
q(t) = 1 9 Steg 


=0 otherwise 


p(t) - q(t) P(f) * Q(f) 


pler ets Pif) * R(f) +f 
Fig.12. Comparison of rectangular and gaussian windows in time and frequency domains 
(a) Infinite sinusoid 
(b) Rectangular weighting function 
(c) Sinusoid weighted with a rectangular function 
(d) Gaussian weighting function . 
(e) Sinusoid weighted with a gaussian function 


used in digital and time compres- When a time window is used in 
sion analyzers since a cosine wave- conjunction with an analyzer filter, 
form or table is already available. the overall filter characteristic is de- 
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Fig.13. Comparison of the spectra of some common time window functions 


termined by the poorer of the two 
with respect to both bandwidth and 
selectivity. (Fig.12). When the filter 
and weighting function have equal 
bandwidth, the combined  band- 
width is about 30% greater (Ref.5). 


Forthe Gauss Impulse Multiplier 
Type 5623 in combination with the 


Analyzer Type 2010, the selectivity 
will virtually always be determined 
by the analyzer (Ref.5), provided its 
bandwidth is greater than that of 
the gaussian window (as it nor- 
mally should be — otherwise the 
analysis takes longer without giving 
an appreciably better resolution). 


This requirement is expressed by 
the following: equation: 


BT, 20,66 (10). 


where 

Tcg= half amplitude length of gauss 
function (indicated on front 
plate of 5623). 


Application to the Analysis of Non-stationary Signals 


It is only when short signal seg- 
ments are to be analyzed that it 
really becomes necessary to use 
time windows. Such is the case 
with non-stationary signals, which 
are segmented into approximately 
stationary sections, and where it Is 
no longer possible to improve selec- 
tivity by increasing the sample 
length. 


Speech analysis (Fig.14) provides 
a good example, where a typical 
sample length is 50Oms. Inserting 
this as T in Equ. (10) it is found 
that the minimum useful bandwidth 
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Amplitude-time function of the English word ’’this’”’ 


200 


Gaussian weighted function used as time window 


300 400 Time ms 


Output from Gauss-multiplier 


171453 


Fig.14. Use of a gaussian window in speech analysis 
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is 13Hz. In any case if a lower 
bandwidth is used (as for example 
happens with 1/3 octave analysis 
for frequencies less than 57Hz2) it 
should be realized that the resolu- 
tion will not be better than 13 Hz. 


Moreover, the minimum valid fre- 


Trigger 


Signal to be analyzed 


After Trigger Recording 


quency is also equal to the mini- 
mum bandwidth. It is in fact worth- 
while to derive the very simple and 
fundamental relationship between 
percentage bandwidth and the num- 
ber of periods in the effective record 
length Te since this gives an im- 
mediate impression of the valid fre- 


0,6 (120 ms) 


quency range for each analysis. (T¢ 
for weighting functions other than a 
rectangle is assumed here to be 
equal to the reciprocal of the effec- 
tive bandwidth listed for the respec- 
tive weighting function, and is ap- 
proximately equal to the half ampli- 
tude length). | 


0,9 (180 ms) 


7502 memory length 200 ms 


Fixed delay 5 ms 


Fixed Delay 5 ms 


Total pulse length 150 ms 


Half amplitude length 50 ms 


Sync. Trigger 


1,2 (240 ms) etc. 


1 Increment , Increment Cc D7 
20 ms 
60 ms = 


Sync. Trigger 


120 ms (etc.) 


Signal sample advances 60 ms 
for each recording on 7502 


(After Trigger recording 0,6 . 
0,9 1,2 ...000 etc) 
Fig.15. 
Thus B woes 
rare 
B__1 
from which oo Tet ; 
b= 11 
or — ( ) 
where 
f = any frequency 
B= minimum bandwidth  deter- 
mined by the weighting func- 


tion 


b= minimum relative bandwidth 
at frequency f. 


n= number of 
quency f in time Te. 


Fixed Delay 5 ms 


periods of fre-. 


Increment , Increment 


Sync. Trigger 


For each recording on 7502, the 
Gauss window is played 3 times 
with Incremental Delays of 0,20 
& 40 ms resp. 


Thus, for b = 1%, n must be at 
least 100 periods etc. 


For example, a 7502 with 10K 
memory includes 2,5K periods of 
the maximum allowable frequency 
(lowpass filter cut off = sampling fre- 
quency/4) and thus includes only 
2,5 periods of a frequency 3 dec- 
ades lower (minimum bandwidth 
40%). If the signal sample only fills 
a fraction of the memory length, 
then the frequency range is reduced 
still further. If it for example corre- 
sponds to 10% of the memory (1 K) 
the frequency range will be reduced 
by one decade. (This will be most 
critical in the case of constant pro- 
portional bandwidth analysis, since 
constant bandwidth gives in itself a 
restricted frequency range). 


a ol eg ee 


Sync. Trigger 


Sync. Trigger 


740327 


Efficient means of moving the time window along the signal 


For this reason it can be an advant- 
age in the analysis of long signals 
to have them recorded on a tape re- 
corder, and either analyze directly 
from a tape loop or (if this results in 
too long an analysis time) play over 
each segment separately to a 7502 
so that it fills as much of the mem- 
ory as necessary. Both ‘After Trig- 
ger Recording” on the 7502 and 
“Incremental Delay’ on the 5623 
can be used to move the time win- 
dow along the signal (Fig.15). 


Independent analyses. will be 
achieved by moving the window a 
step equal to its effective length T_, 
but it can be preferable to reduce 
the step length so that the develop- 
ment of the spectrum can be fol- 
lowed more easily. 
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Analysis of Signal Segments and Transients 


Once a time window has been 
used to select a sample from a con- 
tinuous record, the sample is effec- 
tively an impulse and is analyzed in 
the same way as any other tran- 
sient. Accordingly, this section 
deals with the analysis of transients 
in general. | 


There are two main procedures 
for analyzing transients of which 
the first consists in joining it into a 
loop and analyzing it as a repeated 
Or periodic signal. The other con- 
sists in. playing the transient once 
through the filter for each filter posi- 
tion, and then measuring some par- 
ameter which gives a measure of 
the energy content at that fre- 
quency. 


It is normally simplest and safest 
to employ the first method, using 
RMS detection, since the latter al- 
ways sums the energy content irre- 
spective of phase relationships or fil- 
ter bandwidths. Even though the sig- 
nal has been made aartificially peri- 
odic, this only means that the 
power units measured must be con- 
verted to energy by multiplying 
them by the periodic time T, and 
also that the filter bandwidth should 
be greater than 1/T to acquire a 
continuous spectrum. T can often 
be increased (with blank record) to 
achieve a _ suitably narrow’ band- 
width. : 


The averaging time should be a 
number of times longer than the 
repetition time (to damp out fluctua- 
tions) dependent on how concen- 


trated the energy is in the filter out- 


put signal. It is advisable to try with 
a factor of 3 and then increase as 
necessary. 


There are, however, two limita- 
tions viz. crest factor and dynamic 
range which mean that the method 
cannot always be used. As soon as 
the effective length of the impulse 
falls below about 10% of the loop 
length, there is a chance that one 
or the other factor will not permit a 
valid analysis. In this event it is 
“necessary to employ the. other 
method, which in itself has a.num- 
ber of variations. 


The first is applicable to narrow 
(constant) bandwidth analysis 


where the filter response time TR 
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(~ 1/B) is greater than the impulse 
length. In this case the filter re- 
sponse is virtually independent of 
the impulse shape, and it is possi- 
ble to measure virtually any parame- 
ter (peak, average or RMS value) 
which gives a measure of the en- 
ergy content in the filter response 
pulse. Refs.6 and 7. discuss this 
method, as well as indicating how 
the results may be calibrated. 


Each playback of the impulse 
gives a peak on the recording 
paper, but the correct spectrum can 
be interpreted as the upper enve- 
lope (Fig.16). The filter should not 
be swept more than one bandwidth 
per playback. 


In some cases, one cannot be cer- 


tain that the impulse length is al- 


ways shorter than the filter re- 
sponse time, and it is then once 
more preferable to employ RMS de- 
tection. This is for example the case 
with 1/3 octave real-time analysis 
over a large frequency range (Ana- 


lyzer 3347). It is still possible, how- 
ever, to achieve valid results by us- 
ing an averaging time which is 
several times longer than the im- 
pulse length (or filter response if 
this is longer) but a number of 
times shorter than the repetition 
time, so that the output level has 
time to fall sufficiently between suc- 
cessive playbacks. 


Fig.17 from Ref.2, shows the 
range in which the averaging time 
can be chosen (as a function of the 
filter output pulse length). The area 
where the characteristic is close to 
the “ideal” line can be seen to be 
strongly dependent on the crest fac- 
tor capability of the RMS detector. 
Thus for a crest factor of 5 (e.g. 
3347) the factor Ta/Tg can lie in 
the range 3 to 30. Fig.17 is based 
on rectangular pulses, but should 
be approximately valid for other 
pulses compared on the basis of 
their effective length Te. It is also 
based on the assumption of RC inte- 
gration with B & K detectors. 
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Fig.17. Valid range for the factor Impulse length: Averaging Time (TE/T,a) for RMS detectors 


with different crest factors C 
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For unknown impulses it may be 


difficult to judge Tr, 
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Hat amplitude zone 


Fig.17 can be used to calibrate the 


result, but it would normally be best 
to calibrate using a known signal, 


a rectangular pulse, or tone 


e.g. 
burst. 


Fig.18. Half amplitude length and total length of an impulse 
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